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Abstract 

We present results from a Monte Carlo simulation of non-compact lattice QED in 3 
dimensions on a 16^ lattice in which an explicit anisotropy between x and y hopping 
terms has been introduced into the action. This formulation is inspired by recent 
formulations of anisotropic QED3 as an effective theory of the non-superconducting 
portion of the cuprate phase diagram, with relativistic fermion degrees of freedom 
defined near the nodes of the gap function on the Fermi surface, the anisotropy encap- 
sulating the different Fermi and Gap velocities at the node, and the massless photon 
degrees of freedom reproducing the dynamics of the phase disorder of the supercon- 
ducting order parameter. Using a parameter set corresponding in the isotropic limit 
to broken chiral symmetry (in field theory language) or a spin density wave (in con- 
densed matter physics language), our results show that the renormahsed anisotropy, 
defined in terms of the ratio of correlation lengths of gauge invariant bound states 
in the x and y directions, exceeds the explicit anisotropy n introduced in the lattice 
action, implying in contrast to recent analytic results that anisotropy is a relevant 
deformation of QED3. There also appears to be a chiral symmetry restoring phase 
transition at ~ 4.5, implying that the pseudogap phase persists down to T = in 
the cuprate phase diagram. 

PACS: ll.lO.Kk, ll.15.Ha, 71.27.+a, 74.25.Dw 
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1 Introduction 



The phase diagram of the superconducting cuprate compounds in the (x, T) plane, 
where x denotes the doping, or fraction of holes per Cu02 unit, continues to be the ob- 
ject of much study, both experimental and theoretical. A schematic version is shown 
in Figure HI p. Around x ~ 0.2, so-called optimal doping, there is a superconducting 
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Figure 1: (Color online) Schematic phase diagram for cuprate superconductors 

dSC phase extending to temperatures as high as T ~ 50K. The superconductivity 
is believed to be an essentially two-dimensional phenomenon, being confined to the 
Cu02 planes, and the gap function A(fc) is characterised by c?-wave symmetry, thus 
having two pairs of nodes on the one-dimensional Fermi surface. For x ~ the 
compound is an insulating anti-ferromagnet, and as x increases in this AFM phase 
the order smoothly evolves into a spin-density wave characterised by a wavevector K 
whose magnitude decreases with x. 

In some sense the "normal" phase is more strange. While in the over-doped regime 
the behaviour is that of a normal metal, namely a conventional Fermi liquid, as x 
is decreased unusual non-Fermi liquid behaviour manifests itself via non-standard 
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T-scaling of transport coefficients such as resistivity and themal conductivity. More 
mysterious still is the "pseudogap" behaviour observed in the under-doped region; 
as one moves out of the dSC phase in the direction of increasing T or decreasing 
X, studies of the spectral density distribution function at fixed momentum show the 
quasiparticle pole of the superconductor (correponding to a well-defined excitation 
of energy A above the Fermi energy) diminish in strength, but the magnitude of 
the energy gap |A| remain non-zero even in the non-superconducting region. This 
spectral depletion can persist up to T ~ 150K |2]. 

It should be stressed that while AFM and dSC both have well-defined order pa- 
rameters, and hence are separated from the rest of the phase diagram by solid lines 
in Figure ^ the status of the dashed lines separating the "normal" phase into three 
regions is currently much less clear. Nonetheless, there have been several attempts to 
formulate a theoretical description of the pseudogap region. A particularly interesting 
programme, starting from established properties of the dSC phase, derives an effective 
theory which resembles QED in 2-1-1 dimensions, but having spatial anisotropy in the 
covariant derivatives [SI E] . The starting point is the Bogoliubov - de Gennes model 
for (i-wave quasiparticles in the dSC phase, which in Euclidean metric (corresponding 
to the imaginary time formalism in many body theory) has action 

S = T ^ - %)4(A;,t^„,)c^(A;,u;„) 

- ^ (^A(fc)4(fc, w„)cl^(-fc, -c^„) - A^(fc)c^(fc,u;„)c_^(-fc, -u;„)j , (1) 

where c^,c are creation and annihilation operators for electrons with spin a = ±1, 
and Un = (2n — l)7rT are the allowed Matsubara frequencies. The function is 
the energy of a free quasiparticle, which thus vanishes for k on the Fermi surface, 
and A(/c) is the gap function, which can be thought of as a self-consistent pairing 
field. The requirement of rf-wave symmetry implies that A(fc) = at four special 
node momenta k = iiKi, with K1.K2 = 0. If we choose axes such that Ki \\ x, 
K2 II and write k = Ki + q, then in the vicinity of the "1" nodes it is possible to 
linearise as 

^j: = VFq. + Oiq^) ; Aik) = v^qy + Oiq^) (2) 
and near the "2" nodes as 

^^^ = VFqy + 0{q^) ; Aik) = VAq. + 0{q^), (3) 

where the parameters vp and va are the Fermi and Gap velocities respectively. 
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The next stage is to define a 4-spinor at the node i: 

^^\q, to) = (c+(^, iu), cli-k, -co), c+{k - 2Ki, u), cl{-k + 2i?i, -cu)) . (4) 

The association of different spinor components with different points in /c-space is well- 
known to workers in lattice QCD familiar with the staggered fermion formulation [S]. 
It is now possible to recast the low energy limit of the kinetic term of in relativistic 
garb: 

S = J J dt^i[-fodt + VF-fA + VAl2dy]'^i + 

where /3 = T~^, ^ = \I'^7o, and the 4x4 traceless hermitian Dirac 7- matrices obey 
{""ffijlu} = 2(5^1,. It is important to stress that there is no reason a priori for the 
anisotropy encapsulated in the ratio vp/v^ = k to he negligble in real cuprates; a 
value as high as ~ 10 is reported in j^. In addition, k, increases with doping fraction 

To give the nodal fermions interactions, it is necessary to discuss the reason for 
the loss of superconducting order. The hypothesis is that the gap function can be 
written as A = Ao(fc)e*^*^'^, where Aq is real, and that superconductivity is destroyed 
because the phase field 6 becomes disordered. In two spatial dimensions non-trivial 
phase disorder arises through the accumulation of vortices, that is, point dislocations 
around which § V9.di = 2mi [H]. The dSC pseudogap transition is thus supposedly 
driven by vortex condensation, which preserves |A| ~ Aq but ensures (e*^) = 0. Now, 
it is possible to exploit the gauge symmetry of (Q) to absorb phase fluctuations of A 
into the phases of c, and hence ^ . We would thus seek a theory of phase fluctuations 
for the ^ fields. However, since A represents a doubly-charged Cooper pair field, it is 
impossible to do thus while maintaining single valuedness, since the phase of \E' would 
change by only vr on circling a vortex. The solution proposed by Franz and Tesanovic 
[H], is to partition the vortices of any particular configuration {6} into two groups A & 
B, and then to associate the phase ^a(^) associated with A vortices exclusively with 
the spin- up electrons, and that of the B to spin-down. It can then be shown |31 1^ that 
the relativistic ^ fields of © couple minimally to the vector-valued difference field 
o-ti = \dfj^{9A — Ob), which thus acts as an effective "photon" in the gauge- invariant 
action which results from replacing in (0) by the covariant derivative = d^+ia^. 

It immediately follows from gauge invariance, which implies that the vacuum po- 
larisation tensor correcting {a^{p)au{—p)) has the transverse form {p^d^u — PAtPi')n(p), 
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that the excitations do not receive a mass due to quantum corrections from the \1/ 
fields. Further arguments have been advanced jHl EI to suggest that fiuctuations in 
the field are themselves governed by the action of 2+ld electrodynamics 

Sphot = ^ J '^^^ dt{d^a^ - d^ttf^Y, (6) 

with the coupling g related to the diamagnetic susceptibility x via g ^ x~K in field 
theoretic terms to the dual order parameter $ for vortex condensation via g ~ ($). 

One way of understanding this is that in the absence of magnetic monopoles 
(which in 2+ld are instantons) which can act as a source or sink of fiux, vorticity 
is a topologically-conserved charge. When the ground state is such that ($) 7^ 0, 
the U(l) global symmetry for which vorticity is a Noether charge, i.e. the timelike 
component of a conserved current V^, is spontaneously broken, resulting in a massless 
boson in the spectrum via Goldstone's theorem. However, in 2+ld this Goldstone 
boson is kinematically equivalent to the photon, as can be seen via, eg, the PCAC-like 
relation ^H] 

(OlV^^ll photon, p) oc p^. (7) 

As a result of these considerations, it is natural to unite (0) and © and postulate 
a relativistic field theory, QED3 with Nf = 2 fiavours of nodal fermion as the ap- 
propriate effective action for low energy long wavelength excitations in the pseudogap 
phase - the photons interacting with the \& with effective electric charge g. 

QED3 is an asymptotically-free theory, which means that it becomes more strongly 
interacting as energy scales decrease, ie. as length scales grow. The infra-red be- 
haviour of QED3 has long been a challenge to theory (see JT] for a brief review) ; in 
brief, the issue is whether the chiral symmetry \—>- e^"''^^'^, ^ 1— > ^e^°''^^ of (0) (with 
75 = 7o7i7273) is spontaneously broken by a parity-invariant fermion - anti-fermion 
condensate ('I'\E') 7^ 0. This is believed to happen if the number of flavours Nf is 
smaller than some critical value Nfc, which has been variously estimated as taking 
values in the range ^ | to ~ 5. The consequence is a dynamically generated fermion 
mass E; the determination of the exact value of Nfc, and the dynamically gener- 
ated ratios S/^f^ and (^'\E')/5'^ remain outstanding problems in non-perturbative fleld 
theory. 

To understand this issue in the context of cuprates it is necessary to return to the 
original electron variables c, c^. If Nfc > 2 then chiral symmetry is broken at T = 
in the long- wavelength or continuum limit. This translates into a non- vanishing 
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value for {^^^cr cos {2 Ki.r)cl{r) Ca{r)), which is nothing but the order parameter 
for spin density waves characterised by wavevectors 2Ki This imphes that for 
sufficiently small T there is a direct passage from AFM to dSC without going through 
an intermediate strip of the pseudogap phase, and a corresponding triple point at the 
intersection of AFM, dSC and normal phases for some T > 0. If, on the other hand, 
Nfc < 2, then we expect the pseudogap phase to persist all the way down to T = 0, 
as sketched in Fig. ^ In either case, it may be possible to explain the non-Fermi 
liquid properties in terms of a non-perturbatively large anomalous scaling dimension 
for in the chirally-symmetric phase of QED3 

The above discussion ignores the effects of the anisotropy vf/va = 1^ > 1- This has 
been justified by analytic treatments performed for small departures from isotropy 
in the limit of large- A^j IT^ . where it is shown anisotropy is irrelevant in the 
renormalisation group (RG) sense. Specifically, Ref. J2] uses a Schwinger-Dyson 
approach in the large- A^/ limit to study the behaviour of k under RG fiow, finding 

as Svr^Aj 

where s is the logarithm of the ratio of UV cutoff to physical momentum scales. This 
implies that for weak anisotropy (k ^ 1) Kyen is driven to 1 under RG flow and hence 
irrelevant parameter, ie. '^''^"-^ ^ < 1- Hence, it is argued, predictions 
from isotropic QED3 can be applied directly to cuprates. 

The purpose of the current paper is to examine this claim for arbitrary k, and 
the "physical" case Nf = 2. The theoretical tool we use is lattice simulation of 
so-called non-compact QED, modelling the fermi degrees of freedom as staggered 
lattice fermions. The lattice method is fully non-perturbative, with systematically 
improvable errors due to a flnite UV cutoff (the inverse lattice spacing a~^), and 
a non-zero IR cutoff (the inverse lattice size L~^) of a completely different nature 
to other approaches. The non-compact nature of the model, to be deflned more 
fully in Sec. El below, has the effect of suppressing monopoles (which appear as point 
singularities in the phase fleld a^), thus maintaining the masslessness of the photon 
[HllTn]. Lattice simulations of isotropic QED3 have in the past been applied to the 
issue of Nfc'. such attempts have been hampered by large flnite volume effects in the 
continuum limit due to the massless photon, and to date the results are only able 
to suggest Nfc > 1 CSl- In this study we work away from the continuum limit 
at a coupling g sufficiently strong that we are confident chiral symmetry is broken 
at K = 1, and then systematically increase k. Further details of the lattice model 
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and our simulation are given in Sec. El In Sec. El we present numerical results, and in 
particular present evidence firstly for a restoration of cliiral symmetry at some critical 
Kc, and secondly for the renormalised anisotropy K,ren{i^) > ^5 where Kren is defined 
in terms of certain correlation lengths in differing directions, implying in contrast to 
the analytic results that anisotropy is a relevant deformation of QED3. Implications 
for cuprate superconductivity are discussed in Sec. jH 



2 The Model and Simulation 

The lattice formulation of QED with non-compact gauge fields and staggered lattice 
fermions is described in detail in ref. [ISj. The following is an flavour staggered 
fermion action for QED3 with explicit x — y anisotropy: 

i=l x,x' x,fi<u 

The fermion matrix M^^x' is defined as follows: 

1 ^ 

M^^^> = — ^ ^^,ix) - ^'.x-a^^I'm] + '^^1^'^ (10) 
n=i 

with given by 

= A^r/^(a;) (11) 

and T]^{x) = (— l)^i+---+2:n-i^ where Xi = x, X2 = y and X3 = t, is the Kawomoto-Smit 
phase of the staggered fermion field. The are anisotropy factors, to which we 
assign the following values: = \y = , \t = 1. The purpose of the phase 
factors is to ensure that in the isotropic limit A = 1 the action describes relativistic 
covariant fermions 

If the photon-like degree of freedom 6'^(x) is defined on the link connecting site x 
to site X + fi, then U^^ = exp(ia^a;^) in (fTUI) is the parallel transporter defining the 
gauge interaction with the fermions, and we have a non-compact gauge action given 
by 

Q,,ix) = ^[A;e,ix)-Ate,ix)]. (12) 

The dimensionless parameter (3 is given in terms of the QED coupling constant (ie. 
the "electron charge") via (3 = 1/g'^a, where a is the lattice spacing. It is convenient 
to work wherever possible in units such that a = 1. As discussed above, we use a 
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non-compact formulation of the gauge fields because flux symmetry is not preserved 
in compact U(l) formulations due to instanton formation, which causes the photon 
in such formulations to be massive . 

If we restrict our attention to that portion of the action involving the fermion 
fields, then we see that the introduction of the factors has the effect, at least 
at tree level, of rescaling the lattice spacing in the various directions as = a/ko, 
ay = a/ y/n, at = a. In orthodox lattice QCD similar anisotropics are often introduced 
for technical convenience; eg. spectroscopy of highly excited states such as glueballs 
is considerably more efficient if at -C arc, ay, a ^ EI- ^^^^ ^^^^ ensure Lorentz 
covariance of the continuum limit it is important to check that all terms in the lattice 
action are formulated with the same anisotropy, which results in a fine-tuning problem 
once quantum corrections are introduced. For instance, implementing this programme 
for the action © would require the introduction of separate gauge couphng constants 
Pxt, Pyt, and Pxy, with a non-trivial constraint resulting from the physical requirement 
that eg. the speed of light for photons is the same as that for fermions. In the case 
at hand, though, the plaquette coupling /3 is defined the same in all three planes. 
It is important to stress that in this case the x — y anisotropy is physical, and that 
eg. the resulting ratio a^/ay is an observable to be determined empirically. At 
tree level a^/ay = k; in what follows (see Sec. IH.H.2j) we define this ratio as the 
renormalised anisotropy K^en and estimate it from the spatial decay of a mesonic 
correlation function. Rather than keep track of the various lattice spacings, we prefer 
to think of K as a parameter of the model which can be renormalised through quantum 
corrections. This approach was pioneered in Ref. ^j, where it was shown using large- 
Nf arguments that ^ren < (note that in [12] the equivalent parameter is called A, 
and that our model sets their parameter 5 to 1). 

We set = 1, which yields Nf = 2 fermion species in the continuum limit. An 
algebraic transformation exists relating the single component staggered fields X to 
four-component continuum spinors \l/ fHl, and in particular the chiral condensates are 
related via (xx) = However, we note that in the simulations presented in 

this paper we are working at a strong coupling (/? = 0.2), far from the continuum limit; 
this was done so that we could be reasonably confident that chiral symmetry is broken 
allowing us to examine the effects of anisotropy on the model's phase stucture 
starting from the putative AFM phase. A certain amount of caution is mandated in 
applying our results to the condensed matter-inspired QED3 model ()5l6p . which is 
derived and justified in continuum terms. Further caution is warranted as the fiavour 
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structure of (fTUjl does not entirely capture the theory of jHl Ej; in the condensed 
matter- inspired theory © the second flavour has a vp'^i term in the y direction 
and a f a72 term in the x direction so the two flavours have opposite anisotropics, 
reflecting the fact that there is no physical anisotropy in the original crystal: in our 
model by contrast, following the transformation to \E', 'I' variables the the velocity-7 
matrix structure of the first flavour would be repeated in the second, so that there is 
an overall anisotropy. We expect however that enough similarities between © and 
the cuprate-inspired model persist for us to make reasonable conjectures as to the 
behaviour of the latter system. The point will be discussed further below. 

In order to aquire our results, we utilised a Hybrid Monte Carlo simulation of 
the action with even-odd partitioning on a 16^ lattice aX (3 = 0.2, simulating for 
anisotropics 1 < k < 10 at bare masses of 0.01 < m < 0.05. Further details can 
be found in JH]; here it suffices to note that the algorithm generates representative 
configurations of {9} weighted according to the action © in an exact, that is to say 
unbiased manner. It is in principle possible to perform simulations with a lattice 
action corresponding more closely to the anisotropy structure of but in this case 
simulations would have to be performed with a Hybrid Molecular Dynamics algorithm, 
and results would thus contain a systematic dependence on the timestep size used 
This algorithm would then approximate the "correct" model via a functional 
measure [detM{K)M(K~^)]^; however, away from the continuum limit it remains an 
unresolved issue whether the resulting dynamics is that of a local Lagrangian field 
theory. 

Around 1000 trajectories of mean length 1.0 were generated for each data point, 
and acceptance rates were generally in the region of 70-80% for 0.02 < m < 0.05 and 
60-70% for m = 0.01 (where simulations were less efficient) apart from k = 10.00, 
whose acceptance was over 80%. 

3 Numerical Results 

We performed measurements of the following parameters in our simulation: the mean 
gauge action (f 6^,^) separately in each fiu plane, the chiral condensate (xx)? the lon- 
gitudinal susceptibility xi = d{xx) /dm, and the pion correlator C-,, in each direction, 
from which we extracted the pion mass (t direction) and the effective masses (also 
known as inverse screening lengths) in both space directions x and y. 
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(a) 



(b) 



Figure 2: (Color online) Plaquette action for (a) m = 0.01, and (b) m = 0.05 as 
functions of anisotropy n 



3.1 Plaquettes 

Fig.|21shows average gauge action values for k between 1 and 10 for bare fermion mass 
m = 0.01 (a) and 0.05 (b). This is an important observable in exposing dynamical 
fermion effects: since the plaquette term oc in (jH]) is ^-independent, any anisotropy 
effect must be due to the effects of quantum corrections due to the fermion sector. 

The general upward trend of the x — t and y — t plaquette values with k could 
plausibly be explained as a result of reduced screening of bare charge due to quan- 
tum corrections. In perturbation theory the dominant screening process is known as 
vacuum polarisation — virtual light pairs decrease the effective value of g and 
hence increase the effective jS; this implies that fluctuations (B^) are more strongly 
suppressed by the dynamics of Q. This can be seen by comparing the m = 0.05 and 
m = 0.01 data at k = 1. The increase of (B^) with n would therefore suggest that 
light \E"I' pairs become less important as anisotropy increases. It will prove difficult, 
at first sight, to reconcile this observation with results of Sec. 13.21 

That the value of the average x — t plaquette is consistently less than that in the 
y — t plane can be explained as being due to the overall x — y anisotropy of 0; in 
the condensed matter model (jSJ anisotropic effects should cancel between the two 
flavours. Accordingly, we can interpret the ^ 4% mismatch between (B^j) and (B^j) 
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in Fig. 121 as some measure of the systematic error in our treatment. 

However, the x — y plane plaquettes are markedly different. First, the behaviour 
is non-monotonic - we observe two regimes, one for 1 < k ^ 1.5, and another for 
K ^ 1.5. Within the latter, we have behaviour consistent with the other two planes, 
but for small k the mean plaquette value decreases as k is increased. Secondly, the 
relative splitting between x — y and x — t, y — t is much larger, 0(10%). In this 
case there is no symmetry argument to suggest that this splitting should vanish for 
dynamics based on the condensed matter model ()5I6|) and indeed, the approximately 
quadratic behaviour for small k suggests that some residual anisotropy effect should 
survive even in an x — ?/ symmetrised model, in which effects linear in k might be 
expected to cancel between the two flavours [TTIj . 

3.2 Restoration of Chiral Symmetry 




Figure 3: (Color online) Chiral condensate (xx) versus k. Lines denote fits to the 
equation of state (fT3|) 



11 



The chiral condensate is defined in terms of tlie inverse fermion matrix: 



{xx) 



1 dlxiZ 
V dm 



V 



1 



(trM^^), 



(13) 



and tlie longitudinal susceptibility as: 



Xi = 



d{xx) _ 1 



[(trM"HrM-^) 



{trM-y - {trM-^M'^)]. 



(14) 



dm V 



Anisotropy effects observed within the fermion sector should also be present in the 
symmetrised model although of course with the roles of x and y reversed for the 
second flavour. 

Figures 01 and 0] depict (xx) ^-iid xi fo^' bare masses between 0.01 and 0.05 for 
various k. Note first of all that for k ^ 1 (xx) varies very little as m decreases, and 
certainly extrapolates to a non-zero value as m — 0, implying the spontaneous break- 
ing of chiral symmetry in this limit. This is consistent with the behaviour observed 
at strong couphng in [121 El . Figs. 01 and El both suggest a chiral symmetry restoring 
phase transition as k is increased: a drop in the value of the chiral condensate in Fig- 
ure El (most pronounced for m = 0.01) and, in Figure El a peak in the susceptibility 
which grows more prominent as the mass is decreased. 

We performed a fit of the {xx) data to a hypothetical equation of state of the 
form 



which assumes a second order phase transition at m = 0, k = with conventionally- 
defined critical exponents 6 and Pmag = (^ ~ p)~^ [lEI- We chose to fit values of k 
from 1 to 7, as there are too few points above k = 7.00 to give a good definition of 
the curve. Fixing the values of the exponents p and 5 to 1 and 3 respectively gives us 
a mean field approximation; we find A = 0.0710(2), B = 1.382(6) and Kc = 5.018(7) 
with a x^/d.o.f. of 72. If p and 5 are allowed to vary, we obtain A = 0.103(1), 
B = 2.75(6), = 4.35(2), p = 1.297(9) and S = 3.99(3) with a x^/d.o.f. value of 51. 
The latter fit is plotted as solid lines in Fig. El with the dashed lines denoting the 
equation of state in the chiral limit. Despite the large values of x^, these curves seem 
to describe the data reasonably well; we conjecture that if a phase transition occurs, 
it will do so at ~ 4.5. A finite volume scaling study is needed before the order of 
the phase transition, and hence the validity of (fT3|) can be established unambiguously. 

While Figs. 01 and show clear evidence of a phase transition across which {xx) 
decreases dramatically, it is too early to draw conclusions regarding its precise nature. 



m = A{k- K,){xxY + B{xx)\ 



(15) 
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Figure 4: (Color online) Xi versus k. 

For a second order transition, xi should diverge at the critical anisotropy in the chiral 
limit m 0. However, without a comparison of data from different volumes a first 
order transition cannot be excluded. In the current context the implications are 
profound: a second order transition would lead us to expect {xx) = for m = 0, 
K, > Kc, implying Nfc < 2 in this regime, whereas if the transition is first order it 
remains conceivable that (xx) 7^ is small but non-zero, and hence Nf^ > 2. Caution 
is required because simulations of isotropic QED3 with Nf = 2 cannot exclude a very 
small dimensionless condensate < 10~^ in the continuum limit jTT], invisible 

on the scale of Fig. IHl but nonetheless perfectly consistent with recent analytical 
estimates fK)\ . 

3.3 Pion Correlation Functions and Spectroscopy 

In this subsection we focus on the correlation functions 

C^f^ix^) =^^{x^x{Q)x^x{x)), (16) 
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where the phase e{x) = (—1)^''^''. In the isotropic hmit k = 1 on a symmetric 
lattice all the C-^^ coincide. When chiral symmetry is spontaneously broken it can 
be shown that the correlator is dominated by one of A^^ pseudoscalar approximate 
Goldstone boson poles whose mass oc m. By analogy with particle physics we refer 
to such states as pions; in continuum notation they are interpolated by the operator 
^75^ EI]- Now, in Euclidean quantum field theory for sufficently large separation 
\Xfj_\ the correlator can generally be fitted by the form 

C^^ix^) = A(e-™--"^ + e-'"-''(^--"'^)), (17) 

where is the extent of the lattice in the fi direction. For fi = t the decay parameter 
rriT^t is the pion mass, ie. the excitation energy to create a pion at rest, whereas for 
H = x,y the corresponding quantities are identified as inverse screening lengths. Of 
course, on an isotropic symmetric lattice corresponding to T ^ all three coincide, 
but in our system with explicit x — y anisotropy, m^^, m^^ and m^^ are all distinct. 

3.3.1 Pion correlators in the time direction 

Figure El shows the variation of m.,^t as k is increased. These values were extracted 
from the timeslice pion propagator pfij) via least squares fitting to ()17|1 . Some caveats 
must be offered regarding this data: it is apparent that for the very lightest pions 
the masses were too small for the lattice size (ie. L ^ is not satisfied), meaning 
that we could not use an effective mass plot in order to estimate the ideal fitting 
window for each mass. Instead, we chose the fit window that provided the best 
value and that produced a curve that passed through the error bars of as many of 
the propagator data points as possible. Because of this, a table listing all the data 
points, their values and their fit windows is given below (table HJ. Despite the 
uncertainty this procedure yields in the absolute values of m.„t, the trends found in 
the data are not artefacts of the fit window chosen, coinciding for masses where the 
fit window is more or less stable (m = 0.01) as well as those where it is less so (eg. 
m = 0.05). 

It is clear from Figure El that m^^ for all m increases with k is increased, most 
dramatically for k ^ 4, where the m = 0.01 data show a perceptible kink. It appears 
therefore that in the chiral limit m — > we have two regimes, one where m^t is rel- 
atively insensitive to anisotropy, and one where m^t increases approximately linearly 
with K. It is tempting to identify the boundary between these two regimes with Kc of 
the last section - indeed, non-analytic behaviour across a phase transition should be 
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Figure 5: (Color online) rrint as k increases. 

expected since the pion is a Goldstone mode of the broken phase n < Hc- The linear 
behaviour of m.;^t{K) for all masses in the region 5 ^ k ^ 7 is not as yet understood. 

The behaviour of pions at low energies can be described by the non-linear sigma 
model: 

Snlsm = yJ dHd^.UYid'^U), (18) 

where U (x) = exp(i7r(x) / /^) is a unitary matrix of the chiral group G and the coupling 
known as the the pion decay constant, parametrises the strength of pion self- 
interactions, which become weak in the limit A; — > 0. 

We may calculate for various n by relating and (xx) using the Gell-Mann- 
Oakes-Renner relation [22], 

^If^ = Mxx) (19) 

The results are plotted in Figure IHl They track the (xx) results of Figure El very 
closely; this is perhaps unsuprising, as the chiral condensate was used in our calcula- 
tion, and m,rt varies little for k < n^- 
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Figure 6: (Color online) f-^ as k increases. 
3.3.2 Pion correlators in the x and y directions 

Figure [7| summarises the variation in the effective pion mass in the x and y directions, 
as obtained by fitting Cnx,y data to the form (fTTj) . We were unable to fit for niny 
beyond k = 3.00; the propagator took on a saw-tooth form consistent with the pion 
correlation length being infinite for all practical purposes, ie rriT^y <^ ^y^- Tables |21 
and 01 give fit windows and for each point. 

The data shows m^j, increasing with k, and m^j^ decreasing. This is not unex- 
pected; naively restoring explicit factors of lattice spacing we expect m^^ = M^ra^ = 
A^^MttO where M^r is the expected dimensionful pion mass assuming no physical ef- 
fect as a result of aniostropy. This implies that the mass pole of the propagator is 
shifted by a factor of in the x direction and by in the y direction. Indeed, 
the geometric mean y/ft%i^^frT^ is approximately independent of k, suggesting that 
the results can be explained entirely in terms of equal and opposite anisotropics, ie. 

2 

^y — ^ * 

However, we should consider the possiblity that as a result of dynamical effects the 
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mass(bai"e)=0.01 
□— □ mass(bare)=0.02 
mass(bare)=0.03 
A-A mass(bare)=0.04 
mass(bare)=0.05 




Figure 7: (Color online) Screening masses m^^, (solid), m^ry (dotted), and the geo- 
metric mean ^JnTi^^ffT^ (dot-dashed), versus k. 

physical anisotropy is not simply related to the "bare" anisotropy introduced in (jH)). 
We can then regard the ratio ^Tr^/^Try as a measure of the physical or renormalised 
anisotropy /tren Figure IHlplots the resulting Kren{t^)] in fact Kren is approximately 
described by 



a relation which appears remarkably insensitive to the fermion mass m. 

Eqn. ()2()|1 implies that the physical anisotropy is greater than that in the bare 
theory, in direct contradiction to the analytic prediction of Lee and Herbut jT2] • Our 
results suggest K^en — 1 is relevant. Possible explanations for the discrepancy are firstly 
that we are not necessarily simulating at a small enough anisotropy, or sufficiently 
close to the continuum limit, where the results of jT^j apply, and secondly, as we have 
stressed in Section |2l the model © does not quite reproduce the theory examined 
there. 
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Figure 8: (Color online) Kren versus k in the accessible regime k < 3. 

4 Discussion 

In this paper we have for the first time presented simulation results for a condensed 
matter- inpsired version of lattice non-compact QED3, with a physical number of 
fermion flavours Nj = 2, in which anisotropic fermion hopping in the spatial direction 
has been explicitly introduced. Our main result is that the renormalised anisotropy, 
which we define as the ratio of the pion correlation length in the y direction to that in 
the X, is greater than the bare anisotropy parameter n, and hence that k is a relevant 
parameter in the renormalisation group sense as momentum scales flow towards the 
infra-red. Since the ratio vp/v^ is known to depart from 1 for real compounds, this 
result implies that apparently universal results for, eg. Nfc obtained from QED3 in the 
isotropic limit k = 1 must be treated with caution when applied to superconducting 
cuprates. 

Two caveats should be issued: first, as repeatedly stressed, our model has an 
overall physical x — y anisotropy, whereas anisotropies in the nodal fermion action 
(0) cancel between flavours 1 and 2. We have gone some way towards quantifying 
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this effect with our resuhs in Fig. |2l which show that the extra anisotropy effects 
introduced by our formulation, as manifested by the difference between (0^J and 
(6^j), are significantly smaller than the splitting between these observables and (0^^), 
which must persist in both models. As discussed in Sec. it is possible to formulate 
a lattice model with symmetrised aniostropies to check this issue further, but at the 
cost of using an inexact simulation algorithm. Secondly, since as a Goldstone boson 
the pion is a distinguished particle, it is possible that definitions of K^en in terms 
correlation lengths of other states may yield a different answer. This will be explored 
in future work. 

An interesting and to some extent unexpected result of our study, encapsulated in 
Figs. 121 and m is that there appears to be a chiral symmetry restoring phase transition 
at Kc — 4.5. Strictly speaking, a finite volume scaling study on a range of lattice 
volumes will be needed to elucidate the nature of the phase transition, but if it proves 
to be second order then it is difficult to avoid the conclusion that = for large 

anisotropies even for Nf = 2, implying that Nfc is a decreasing function of n, and that 
therefore the pseudogap phase persists down to T = in cuprates. Physically, the 
phase fluctuations hypothetically responsible for the destruction of superconducting 
order must then arise as a result of quantum, as opposed to thermal, effects. 

We can estimate the range of values of k, for which our results might in fact be 
physically relevant. The empirical equation for the boundary of the dSC region [21] 

^ = l-82.6b-0.16f, (21) 

c 

used for the cuprate YBCO, where p is the hole concentration, gives us a value of 
p ~ 0.05 for Tc = in the underdoped region. Sutherland et al. have measured k 
for four values of the doping of this cuprate; from Figure 4 of [7j one can extrapolate 
by eye that 6 ^ k ^ 8 at the onset of the superconducting phase. This does seem to 
suggest that the QED model predicts the occurence of a phase transition at a value 
of K within the region of its validity; were Kc > 6 the phase transition would occur 
after the onset of superconductivity, which would be unphysical. 

Finally, it is interesting to speculate on the nature of the chirally symmetric high-K 
phase. Franz et al P] have argued that in the chirally symmetric phase the fermion 
propagator (\&(0)^^(x)) receives a large anomalous scaling dimension from quantum 
corrections, which are calculable as a power series in NJ^. We find it difficult to 
reconcile these ideas with the plaquette data of Fig. |21 which show plaquette fluctu- 
ations increasing with /t, whereas massless fermions would be expected to suppress 
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such fluctuations through screening. A plausible alternative is that chiral symmetry 
restoration in the model is itself driven by fluctuations of the phase of as the 
system dynamics becomes more and more two dimensional with increasing k, but 
that fermion mass generation, which depends on (l^^^l), remains insensitive to k. 
Similar effects have been observed in model simulations in which phase fluctuations 
are thermally driven [2^] • In future work we intend to explore this issue further with 
measurements of the gauge-fixed fermion propagator. 
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Table 1: Pion mass 171^^1 fitting data 



22 



m 


K 




d.o.f 
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0.01 
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0.289(1) 
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0.742 
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0.975 
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1.735 


3-13 


0.02 
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0.298(1) 


1.256 


1-15 
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1-15 
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Table 2: Effective pion mass m^x in the x direction 
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Table 3: Effective pion mass m^ry in the y direction 
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